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Abstract 

A theory for predicting the noise field of supersonic 'propellers with realistic blade 
geometries is presented. The theory, which utilizes a large-blade- count approximation, 
provides an efficient formula for predicting the radiation of sound from all three sources 
of propeller noise. Comparisons with a full numerical integration indicate that the levels 
predicted by this formula are quite accurate. Calculations also show that, for high speed 
propellers, the noise radiated by the Lighthill quadrupole source is rather substantial when 
compared with the noise radiated by the blade thickness and loading sources. Results from 
a preliminary application of the theory indicate that the peak noise level generated by a 
supersonic propeller initially increases with increasing tip helical Mach number, but it 
eventually reaches a plateau and does not increase further. The predicted trend shows 
qualitative agreement with the experimental observations. 


Introduction 

Over the last dozen years, the renewed interest in supersonic propeller (propfan) 
technology on the one hand, and the drive to meet stringent community and cabin noise 
regulations on the other, have spurred a great deal of theoretical and experimental research 
activity in the area of high speed rotor noise. 

On the theoretical side in particular significant strides have been made towards the 
development of accurate theories for predicting supersonic propeller noise. The foundation 
for almost all of these theories has been the Ffowcs Williams and Hawkings 1 (FW-H) 
equation, which provides a theoretical basis for calculating the radiation of sound from the 
propeller blade surfaces (i.e., thickness and loading noise) as well as the sound generated 
by the flow field surrounding the propeller disc (i.e., Lighthill quadrupole noise).. The 
noise generated by thickness and loading sources has received a great deal of attention as 
exemplified by the theories of Hanson 2 and Farassat. 3 In contrast, most of the research 
in the area of propeller quadrupole noise has been of an exploratory nature (see, for 
example, Hanson and Fink 4 ) even though conclusions from these efforts have suggested 
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that quadrupole radiation could be potentially significant for supersonic rotors. Recent 
work by Peake and Crighton, 5 however, has reinforced the notion that quadrupole radiation 
must be accounted for in predicting the noise generated by high speed propellers. 

FYom a computational point of view, prediction of propeller noise involves evaluation 
of the multiple integrals appearing in FW-H equation. This task can be accomplished 
by various methods but, in general, an accurate and detailed prediction requires the use 
of numerical integration. That, in turn, necessitates the utilization of a large number of 
source points to produce a realistic representation of the propeller geometry as well as the 
flow field surrounding it. In practice, the number of thickness and loading sources necessary 
for such a representation is sufficiently small so as to make the computation of radiation 
by these two types of sources quite feasible. In contrast, a substantially larger number of 
quadrupole sources needed for a faithful representation of flow field around a propeller disc 
makes quadrupole noise computations considerably less viable. While the application of 
various simplification schemes like the far-field approximation has yielded some measure of 
success, a realistic calculation has, for the most part, remained computationally expensive. 

An appealing method for circumventing this problem, first suggested by Hawkings and 
Lawson 6 in the context of a frequency-domain analysis, is the use of the large-blade-count 
approximation to calculate the source distribution integrals appearing in FW-H equation 
asymptotically. For modem, many-bladed propfan designs this approach offers a very ef- 
fective alternative to direct numerical integration. The effectiveness of this approach was 
demonstrated by Crighton and Parry, 7 who developed an asymptotic theory for propeller 
thickness and loading noise radiation. In fact, this theory was also used by the authors in 
Ref. 5 to assess the importance of propeller quadrupole noise. To date, however, the pub- 
lished results of this theory have been limited to propellers with fairly simple geometries. 
Therefore, as it stands, there is a gap between the theories which utilize full numerical 
integration and can accommodate realistic blade geometries but are rather unwieldy for 
calculating noise radiation from quadrupole sources, and the asymptotic theory of Ref. 5, 
which is quite efficient for quadrupole noise calculations but, so far, is restricted to simple 
geometries. 

The theory presented in this paper bridges that gap. It utilizes a large-blade-count 
approximation to evaluate integrals in FW-H equation. But, unlike Refs. 5 and 7, the 
asymptotics is applied to the radiation efficiency integrals instead of source distribution 
integrals. This approach affords the flexibility of using realistic blade geometries while 
keeping the computer time requirements small. The asymptotic analysis for thickness and 
loading noise sources was presented in Envia 8 and shown to provide a very good agreement 
with direct numerical integration. The theory is now extended to include quadrupole noise 
and, therefore, for the sake of completeness the full analysis will be detailed here. Compar- 
isons with a full numerical integration will then be shown and a preliminary application 
of the theory to a practical problem of interest will be presented. 

Analysis 

The starting point for the analysis is Goldstein’s 9 version of FW-H equation for the 
acoustic pressure p(x, t) written in the propeller-fixed (wind tunnel) coordinate system: 
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where v n is the normal component of the blade surface velocity, f is the amplitude of 
the aerodynamic loading on the blade and Tjk is the Lighthill stress tensor, p, p and Uj 
are the fluid density, pressure and velocity with the subscript “o” denoting their ambient 
values. The velocities Uj and v„ are given relative to a medium-fixed coordinate system 
even though Eq. (la) is expressed in a propeller-fixed coordinate system. S(t ) and V(r) 
represent the propeller blade surfaces and volume surrounding the blades, respectively, 
and rij is the component of the outward unit normal to S. It should be noted that 
the convention for the direction of unit normal is opposite that of Goldstein s. y and x 
represent the source and observer coordinates, respectively, t is the source time and t is 
the observer (i.e., retarded) time. U 0j denotes the medium convection velocity which is 
equal to, and negative of, the propeller forward flight velocity. The three terms in Eq. 
(la) represent contributions from the thickness, loading, and Lighthill stress (quadrupole) 
sources of propeller noise, respectively. 

G = G(x, </ y, r) in Eq. (la) denotes the free-space, moving-medium Green’s function 
which can be shown to have the following form: 


G = 7— —5 8 (t — t — g c R/Co) (2a) 

4ttkR 

Sc(r) = 4(« - k(t) = (M 0 2 „ + ft)' 12 (2b) 

P 0 

M 0r (t ) = M 0j ej, Po = (1 - Mq. ) 1/2 (sum over j) (2c) 

Mo,. = e i ( r) = ^( r ) = l x - y( T )l- ( 2if ) 


Here Mq r is the medium convection Mach number in the radiation direction, R is the dis- 
tance between observer and source, ej is the component of the unit vector along radiation 
direction, and C 0 is the medium ambient speed of sound. Parameters g c and k represent 
the effects of medium convection on the retarded time and spherical spreading rate, re- 
spectively. The explicit dependence of various parameters on source time r is indicated 
where necessary. 
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In general, because propeller may be operating at an angle of attack to the oncoming 
flow, the 1-axis is chosen to coincide with the propeller shaft (see Fig. 1) in order to simplify 
the description of the source motion. As a result, motion of the sources is confined to the 
transverse planes described by the 2— and 3— axes. Therefore, it is readily seen that 


Mqj = Mo cos a, Mq 2 = 0, Mo a = Mq sin a 


(3) 


where a is the propeller angle of attack with respect to oncoming flow as depicted in Fig. 
( 1 ). 


Provided that blade geometry, aerodynamic loading, and flow field about the propeller 
are known in sufficient detail, acoustic pressure p(x, t) may be computed directly from Eq. 
(la) for any observer location. This, of course, is the so-called time domain approach 
which entails solving transcendental equations for source time r. Alternatively, one could 
use the frequency domain approach, which involves expanding p(x, t) in terms of its Fourier 
harmonic components p,(x), i.e., 
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/= — oo 

with the individual harmonic components given by: 
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in which ft is the propeller angular speed. P 7 j, pc l and pq ( denote contributions to the 
harmonic component p ; from the thickness, loading, and quadrupole sources, respectively. 
In writing Eqs. (4c-4e) the derivatives of G with respect to source coordinates from Eq. 
(la) have been replaced by its derivatives with respect to observer coordinates through use 
of the relations 
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Using the usual phase relationship arguments it can be shown that, given B identical 
blades, only harmonic components for which l = mB (where m is an integer) contribute 
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to the infinite sum in Eq. (4a). This contribution is simply B times the contribution of 
a single blade. Therefore, from here on B will explicitly appear in expressions for the 
harmonic components p Tmfl , Pc mB , and p QmB . Correspondingly, 5 and V, respectively, 
will denote surface of, and volume around, a single blade. 

To further simplify Eqs. (4c-4e) the spatial derivatives d/dxj and d /dxjdx k of G 
may be rewritten in terms of the temporal derivative d/dt through the use of the chain 

rule: 
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where the auxiliary derivative operator Vj{t) has been introduced for notational brevity. 
The temporal derivatives can now be removed using integration by parts, and the integr 
over t, which involves a delta function, can then be easily evaluated to yield 
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where and represent the expressions for the strengths of thickness, loading 

and quadripole sources, respectively. Note that the Q^sfor each source type are organized 
according to which power of radiation distance R they multiply. Note also that the limits 
on the integration over r have been changed to reflect the fact that an interval of size 
2ir/fl in t is mapped exactly into an interval of size 2r/Sl in t. The expressions for source 

strengths are given by 


Q Q = — — imBQ pov n 1 Mo^ej — g c ^0j) 

1 k 1 K 

Qf = - 1- p 0 v n M 0j (Pfe + M 0R M 0i ) 


(8a) 

( 86 ) 


5 


(8c) 


Qc^ = ~c^2 imB f n i ( e J - 9cM 0j ) 



Qf 

II 

n j (0o e j + Mo r Moj ) 



(W) 


Qa 

1 

~ k'+ 2 

T jk A < ! ) 6j k + A{ t) e j e k +A i 3 e) M 0i 

M 0k 






+ A^ (cjMo* +M 0i e*)j t = 

1,2,3 


(8c) 

4° = o, 

4° 

= 1, 

A {1) - a 2 
^3 9 c 5 

A 1} 

= -9c 


A^ = — K 2 , 

- 4” 

= 3ft 2 , 

— —(1 + 

4 2) 

= 2M 0r - 

Po 9c 

* 

<Np 

1 

II 

e? 

■w «1 

_4 ( 3) 

= 3(ftV 

- MS,), 4 3) = (3 MS, - -c 2 ), 

4° 

3 

s 

<£p 

CO 

II 

. 


m 


where sums over the indices j and k are assumed. It should be noted that, in general, 
normal surface velocity v n , loading amplitude /, and the Lighthill stress tensor Tii depend 
on source time r. 

Once source strength distributions Q^} are estimated, what remains is to carry out 
the integrations in Eqs. (7a-7c). Ordinarily, the order of integration in these equations is 
reversed and the integral over r (which represents the radiation efficiency of the acoustic 
source) is computed first. Since for general geometries and source strengths this integral 
is not tractable analytically, it is usually computed approximately for near- or far-field 
observer locations for which the integrand may be significantly simplified. The result is 
given in terms of the appropriate Bessel functions. The remaining surface (or volume) 
integral is then carried out using a quadrature scheme. Alternatively, as was mentioned 
earlier, surface (or volume) integral may be calculated asymptotically in the manner sug- 
gested in Ref. 6. In this approach, which is particularly useful when the blade count B (or 
more appropriately the harmonic index mB) is large, the method of stationary phase is 
utilized to show that, asymptotically, most of the radiation from blade surface (or volume 
surrounding it) comes from the neighborhood of special points called the stationary phase 
points. As was pointed out in introduction, this approach forms the basis of theory pre- 
sented in Refs. 5 and 7. The utility of this approach hinges on determining, analytically, 
the stationary phase points - a difficult task when dealing with blades having complicated 
geometries. 

This difficulty, however, can be circumvented by applying the large-blade-count ap- 
proximation to integral over r instead. In other words, the idea is to find the radiation 
efficiency of each source asymptotically. This is accomplished by evaluating the r integral 
using a modification to the standard steepest descent and saddle point methods. The 
advantage of this approach is that it is applicable to general geometries. 

In this paper the derivation for the zero angle of attack case ( a = 0) will be outlined. 
The extension to nonzero angle of attack follows the analysis described in Ref. 8. For the 
sake of brevity the final formula will be given in terms of a generic source strength Q and 
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a generic radiation distance “factor” H. The resulting expression is therefore applicable 
to any Q^) and R* dependence. 

For practical considerations it is more convenient to preserve the order of integration 
in Eqs (7a-7c) and carry out surface (or the volume) integral first. To do so, we begin 
by dividing blade surface (or volume surrounding it) into a number of small surface (or 
volume) elements. Let that number be N s , where the subscript 3 indicates a particular 
source element under consideration. If the typical elemental size is sufficiently small, 
integrand may be assumed to be constant over the extent of element and thus surface 
(or volume) integral may be approximated by some appropriate “mean value” of integrand 
times the elemental surface area (or volume). If the mean value is chosen to be the value of 
integrand at the geometric center of element, designated y s , the error in this approximation 
is 0(L 2 ) for a given N a where L = |y - y 8 |n, a x- Of course, error could be made arbitrarily 
small by choosing a large enough N s , which would reduce the effective size of each element. 

Therefore, the generic harmonic amplitude V SmB (x) may be written as: 
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where (x\,X 2 ,x$) and (yi,,y 2 ,,y 3 .) have been replaced by their cylindrical polar counter- 
parts (xi,r, <f>) and (yi,,r ai <f>g), respectively. Furthermore, the integration variable r is 
now replaced by a new variable 9. M t i p is the tip rotational Mach number (i.e., R t i P (l/C 0 ) 
and c, is a convective phase factor representing the collection of phase terms which do 
not depend on r. Q a (9) and 'R a (8) are the generic source strength and radiation distance 
factor, and A a is elemental area (or volume) for the sth element. In order to simplify 
the notation, we suppress the explicit dependence of various variables on y s , the centroid 
of sth element. $,(#) is the canonical phase function for a propeller operating at zero 
angle of attack, with a 3 and b a representing a combination of geometric, convective and 
kinematic factors. Note that in deriving $*(#), no assumption regarding the geometry of 
propeller blade or the location of observer has been made. For this reason, the subsequent 
results apply to arbitrary geometries and observer locations. In the above expressions both 
observer and source spatial coordinates are nondimensionalized by the propeller tip radius 
Rtip. (Note that M 0l = Mq for zero angle of attack case.) 
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Now, for a typical propfan B — 8. Thus even for the blade-passing-frequency (BPF) 
tone (i.e., m = 1) the integrand in Eq. (9b) is highly oscillatory. Nallasamy et al., i0 who 
considered only the thickness and loading noise contributions, employed a numerical inte- 
gration scheme to calculate I 9 . Unfortunately, accurate computation of higher harmonics 
requires an ever increasing quadrature resolution to capture the oscillatory nature of the 
integrand and that can increase the computational cost substantially. The asymptotic 
evaluation of I 9 for large mB, however, can be quite cost-effective. 

In order to carry out the asymptotic evaluation, variable 9 must be allowed to be 
complex and the integration path (i.e., [0, 27r]) be deformed into an appropriate contour. 
Replacing 0 with the complex variable v — 9 + ia leads to a phase function $ a (^) which 
is now also complex. It is fairly straightforward to find the saddle points of $ 9 (v) (i.e., 
$ a (i/) = 0) and the appropriate steepest descent contours. The easiest way to find the 
saddle points is to define an auxiliary variable £ = cosi/ and rewrite $ a (i/) = 0 in terms of 


M,/2) 2 i 2 - b.t + 1 - (a s b 9 / 2) 2 = 0 (10) 

where a 9 and b 9 were defined in Eq. (9d). Eq. (10) is clearly quadratic in £ and may readily 
be solved. Thus for given observer and source locations $(i/) has, in general, two simple 
saddle points in the interval [0, 2n] which have different forms depending on whether the 
component of the source relative Mach number in the direction of observer, A/ r , is subsonic 
or supersonic. The two are a complex conjugate pair if M r < 1 and are real if M r > 1. 
When M r = 1 (i.e., the “sonic condition”) these two saddle points merge to give rise to a 
single, second order saddle point. For a subsonic source, only one of the two saddle points 
lies on the appropriate steepest descent path. Hence only that saddle point contributes to 
the integral. It should be noted that the integrand also has an infinite number of branch 
points coinciding with those of $ a (i/), of which only four lie in the region of interest. A 
judicious choice of branch cuts guarantees that the contributions to the contour integral 
along branch cuts in the I'-plane is exponentially small compared with the contributions 
from the neighborhood of saddle points. The location of the appropriate saddle points, 
branch cuts, and steepest descent contours for typical source locations are shown in Fig. 
2. Note that the choice of steepest descent path depends on whether M r is less than or 
greater than unity. Also shown in this figure are the auxiliary descent contours needed 
for deforming the original contour [0,27r] into steepest descent contour(s). Due to the 
periodicity of integrand in 9, the contributions from these contours cancel each other out 
exactly. 

For a given observer location and operating condition, the asymptotic structure of I np 
for a source with M r = 1 is different from that for a source with M r ^ 1. In fact it, is fairly 
easy to show that the asymptotic expansion is given in terms of the inverse i/3-powers of 
parameter mB for a “sonic” source and in terms of the inverse 1 / 2 -powers of mB for a “non- 
sonic” one. Since there is no simple way of constructing a composite expansion from these 
two expansions to allow for a smooth transition through the sonic condition, they axe not 
very convenient to use. However, this difficulty may be avoided altogether by developing 
a uniform asymptotic expansion using a theorem derived by Chester et al. 11 The details 
of the methodology may be found in Bleistein and Handelsman. 12 The basic idea is to 
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map, conformally (say, v -> (), phase function $ a (u) into a much simpler function (i.e., * 
cubic polynomial) which exhibits the relevant features of the original phase function. The 
region of interest in complex i/-plane (i.e., the region containing saddle points and steepest 
descent contours) is correspondingly mapped into a region in the complex (-plane. The 
standard steepest descent and saddle point methods are then applied to the integral in the 
(-plane. The key definitions and parameters, along with the final result, axe summarized 

below. 

The cubic is given by: 



(11a) 

= 5 [*.(>'?) + *.(01 

(111) 

il = I [M-'?') - *.(Ol 

(He) 


where i/+ and u~ denote the locations of saddle points of $ 8 in complex i^-plane and 
and 7, are parameters defining the conformal map. Note that 7, as given above can take 
on three possible values or branches. The theorem in Ref. 11 guarantees that one of the 
branches defines the desired conformal map. For integral I a the branch chosen is the one 
for which 7* is purely real. With these parameters determined, the map can be constructed 
and the uniform asymptotic expansion carried out. In principle, expansion of integral /, 
could be developed to an arbitrary order in the parameter mB. However, it turns out that 
for most applications the first term provides a very reasonable approximation even at BPF 
tone (i.e., m = 1 ). Therefore, in the subsequent analysis only the first term is considered. 

After a fair amount of algebra the final result can be written as the following formula: 


Is 


~ 2tt 


, Ai({mB)>f 3 7?) , J Aj ((mB)^S?) } 

So (m 5)>/ 3 + 41 (mB ) 2 / 3 j 


d ao 


r 0 (7*) + Fo(— 7 *) 




r 0 (7«) - r 0 (-7«) 

27- 


r rs\ - — 

r ° (C) " K a (v«)) dC 


du _ 7? ~C 2 

d( *>«)) 


(12a) 
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where Ai and A[ axe the Airy function and its derivative, respectively, and d ao and d ax are 
coefficients in the asymptotic expansion. It is worth mentioning that ±7* corresponds to 
the locations of saddle points vf in (-plane. For any observer location, Airy function and 
its derivative provide a smooth transition from the region of blade for which M r is less than 
unity to the region of blade for which M r is greater than unity. For 7« = 0 (i.e., a sonic 
source) A, and A\ are 0 ( 1 ) and consequently I a is proportional to inverse 1 /3-powers of 
mB. For 7, ^ 0 (i.e., a “non-sonic” source) and large mB Aj is 0 ((mB)" 1/s ) and A; 
is 0 ((mB)‘ /s ). Consequently, I a is proportional to inverse 1/2-power of mB as expected. 
After substituting for I, in Eq. ( 9 a) from Eq. (12a) and adding the contributions from all 
the N a elements, Fourier harmonic component 'P m a(x) can be calculated. 
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Results and Discussion 


In order to assess the accuracy of the asymptotic formula given by Eq. (12a), we next 
consider its predictions for a test problem. The test case is that of a propfan operating 
at cruise conditions. The necessary aerodynamic input to the acoustic model is computed 
using an Euler CFD code developed by Adamczyk (see Celestina et al 13 ) which provides 
both the loading distribution on propfan blades and flow field around the propfan. The 
CFD computations are for the flow conditions which roughly correspond to the wind tunnel 
conditions in a series of acoustic measurements carried out at NASA Lewis Research Center 
by Dittmar and Stang 14 for a scale model SR-7A propfan operating at simulated cruise 
conditions. All of the acoustic calculations presented in this paper are for BPF tone noise. 

As was previously mentioned, the asymptotic formula is applicable to all three types 
of sources represented in Eqs. (7a-7c). The application of formula to thickness and loading 
noise sources (i.e., surface sources) is quite straightforward since the spatial extent of their 
distribution is clearly defined. However, the “effective” spatial extent of the distribution of 
quadrupoles (i.e., volume sources) is not known a priori. To determine this effective spatial 
extent a series of calculations employing the asymptotic formula for the quadrupole sound 
pressure level (SPL) was carried out using a series of progressively larger computational 
volumes surrounding the propfan blade. 

The results of this study are shown in Fig. 3 where the predicted sideline directivities 
of quadrupole SPL (at BPF) for different computational volumes are plotted for near-field 
observer positions located at a distance of 1.6J?(, P from the axis of propfan, and extending 
from 1.0J?t,'p forward to 1.0i?t,'p aft of the plane of rotation. The medium convection (free 
stream) Mach number Mq is 0.8 and the advance ratio of propfan is 3.06, which together 
result in a (relative) tip helical Mach number (M tiPh = (MZ+M? ir ) 1/3 ) of 1.15. As a matter 
of convenience the computational volumes were chosen to coincide with the subsets of 
the CFD computational grid. In all, four volumes were considered. The passage volume, 
denoted VO, extends axially from the blade’s leading edge to its trailing edge and extends 
radially from hub to tip. Volume VI extends axially to about 0.6f2(,-p, volume V2 to 2.0 Rti p , 
and volume V3 to 4.0 R t i p in both the forward and aft directions. The latter three volumes 
also extend a distance of 0.5i2t, p beyond the blade tip in the radial direction. All four 
volumes also span the blade passage halfway to the adjacent blades in the circumferential 
direction on each side. For each volume the contributions to noise field from Lighthill 
quadrupole sources contained within the volume were computed and plotted. 

Two important conclusions drawn from the plot in Fig. 3 are that (i) the predicted 
peak quadrupole SPL occurs in the vicinity of plane of rotation and (ii) that as far as the 
peak noise levels are concerned, volume VI is optimally suited for acoustic calculations 
since the larger volumes produce levels which are generally no more than 1 dB different from 
those for volume VI (the results for volumes V2 and V3 are virtually indistinguishable). 
The erratic behavior of predicted SPL for the baseline volume VI at the axial observer 
locations aft of the 0.6Ru p position is a result of the choice of computational volume. 
Because, these observer locations lie near the boundary between VI and V2, in the baseline 
case they only receive contributions from sources located upstream of that boundary (i.e., 
those inside VI only). Clearly, as computational volume is extended and contributions 
from sources downstream of that boundary (i.e., those inside V2 but not inside VI) are 
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also taken into account, more realistic levels are produced and the “kink in the SPL 
disappears. Furthermore, in order to verify that the kink is not an artifact of asymptotic 
calculations, a full numerical integration of radiation efficiency integral for volume VI was 
performed at a few discrete locations and the results plotted as solid symbols in Fig 3. The 
comparison between asymptotic and numerical calculations is strikingly good and serves 
to demonstrate the accuracy of asymptotic calculations. 

To further emphasize this accuracy, a comparison between the asymptotically and 
numerically computed sideline directivity of total SPL (thickness + loading + quadruple) 
is presented in Fig. 4. The operating conditions and observer locations are the same as 
those in Fig. 3. The quadrupole contribution from volume VI was considered in this 
comparison. The numerical integration was carried out for only five observer locations to 
limit the computational time. Again, the agreement is very good with deviations of less 
than 1 dB in plane of rotation between asymptotic and numerical results. For comparison, 
the corresponding sideline directivity for the noise generated by surface sources is also 
shown. The agreement is quite good with a maximum deviation of about 2 dB occurring 
at about 0.5 Ru p downstream of plane of rotation. A detailed examination of the individual 
source type contributions in this region revealed that on the one hand the discrepancy is 
partly due to deviations (of around 2 dB) between the asymptotically and numerically 
calculated loading noise predictions, and on the other to the strong sensitivity of the sum 
of thickness and loading noise SPLs to their individual contributions. It is also interesting 
to note that the noticeably higher total SPLs as compared with those for surface sources 
indicate that quadrupole noise levels are rather substantial for high speed propellers and 
must be taken into account. 

Finally, to demonstrate the utility of the method, a series of acoustic calculations was 
performed to investigate the sensitivity of the maximum total SPL generated by propfan 
as a function of its tip helical Mach number M tiph . The increase in M tiph was achieved by 
holding the advance ratio constant (at 3.06) and increasing the free stream Mach number 
Mq. The range of tip helical Mach numbers studied was between 0.9 and 1.4 with 11 
different Mach numbers considered within this range. The aerodynamic input in each 
case was computed using the Euler code mentioned earlier. As before, the computational 
volume chosen for acoustic calculations was VI. Twelve observer positions distributed 
on both sides of plane of rotation were considered. These positions correspond to the 
microphone locations in Ref. 14. In each case the asymptotically-predicted maximum SPL 
was noted amongst observer positions considered. In all cases this maximum occurred in 
the neighborhood of plane of rotation. 

The peak SPLs are plotted versus M t i Ph in Fig. 5. To see the trend, a quadratic 
curve is fitted to the predicted values (solid circles). The predicted trend indicates that 
maximum SPL does not increase monotonically with increasing tip helical Mach number, 
but rather that it levels off beyond M tiph of about 1.25. This behavior is consistent with 
the trends observed in the measurements reported in Ref. 14. As seen from the inset in 
Fig. 5, measured peak SPLs generally reach a plateau past a tip helical Mach number of 
about 1.25. Different symbols in the inset represent peak SPLs corresponding to different 
propfan blade setting angles while the lines are quadratic curve fits through the points. 
The advance ratio in all three cases is 3.06. In this paper no attempt was made to carry 
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out a detailed comparison with the experimental data due to the ambiguities which exist in 
matching the experimental operating conditions with their CFD counterparts. For the sake 
of comparison, however, the maximum SPLs as predicted from only the surface sources 
are also plotted (open circles) in Fig. 5. It is interesting to note that, in contrast with 
total noise, predicted peak noise arising from surface sources, tends to rise monotonically 
with the increasing tip helical Mach number. An examination of theoretical predictions 
shows that, the difference in behavior between peak SPLs from surface and total sources 
may be traced to the mutual phase relationship between surface and volume sources. For 
M t i Ph < 1.25, surface and volume sources radiate essentially in phase. Therefore, they 
contribute additively to total SPL. For M t i Ph > 1.25, however, surface and volume sources 
radiate out of phase. This leads to mutual cancellation between the two types of sources 
rind, as a result, in reduced peak SPLs. 

It is worth mentioning that the CPU times required for these asymptotic noise calcu- 
lations were generally an order of magnitude smaller than those for full numerical integra- 
tion. The CPU time savings would be even more significant for higher harmonics of BPF 
where numerical integration requires a progressively finer quadrature step size. In contrast, 
asymptotic calculations would incur a nominal increase in CPU time over that required 
for the BPF tone calculations, since saddle point locations are independent of the mode 
number. Furthermore, due to the very nature of the asymptotic expansion, predictions 
for the higher harmonic will be even more accurate when compared with numerical results 
(see Ref. 8). 

Concluding Remarks 

A large-blade-count asymptotic theory which allows for accurate and efficient calcula- 
tion of noise field of high speed propellers was presented. The theory does not rely on the 
simplifying assumptions usually employed in other propeller noise analyses. A closed form 
expression involving Airy function and its derivative gives a uniform representation of the 
noise field of a source regardless of whether its component of the relative Mach number in 
the direction of the observer is less or greater than unity. The levels of all three types of 
propeller noise sources (i.e., thickness, loading, and Lighthill quadrupoles) were computed 
for a realistic propeller geometry and operating condition at a very reasonable computa- 
tional cost. The results indicate that near-field quadrupole radiation is rather substantial 
in vicinity of plane of rotation and can significantly influence total sound pressure levels 
radiated by a supersonic propeller. Furthermore, inclusion of quadrupole contribution in 
predicted levels seems to generate trends which are in qualitative agreement with exper- 
imental measurements of sensitivity of maximum SPL to increases in tip helical Mach 
number. The preliminary results shown here indicate that the present method provides a 
useful theoretical tool in propeller noise research. 
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Fig. (1) Definitions of Geometry and Coordinate System. 
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Fig. (2) Saddle Points and Steepest Descent Contours in the 
Complex i/-Plane. 
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Fig. (3) Sideline Directivity of the Predicted Quadrupole SPL 
as a Function of the Computational Volume Chosen. 
(Asymptotic Calculations) 
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Fig. (4) Sideline Directivity of Thickness + Loading + 
Quadrupole and Thickness + Loading SPLs. 
Comparison of Asymptotic and Numerical Results. 
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Fig. (5) Predicted Variations of Peak Thickness + Loading + 
Quadrupole and Thickness + Loading SPLs with Tip 
Helical Mach Number. Inset: Experimental Data. 
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